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ABSTRACT 

We present high-resolution, //-band, imaging observations, collected with 
Suhani/HiCIAO, of the scattered light from the transitional disk around SAC) 
206462 (lit) 135344B), Although previous sub-nun imagery suggested the ex- 
istence of the dust-depleted cavity at r < 46 AU, our observations reveal the 
presence of scattered light components as close as (f!2 28 AU) from the 

star* Moreover, we have discovered two small-scale spiral structures lying within 
C//5 70 AU)* We present models for the spiral structures using the spiral den- 

sity wave theory, and derive a disk aspect ratio of h ~ U.L wliich is consistent 
with previous sub-mm observations* This model can potentially give estimates of 
the temperature and rotation profiles of the disk based on dynamical processes, 
independently from sub-mm observations. It also predicts the evolution of the 
spiral structures, which can be observable on timescales of 10-20 years, providing 
conclusive tests of the model. While we cannot uniquely identify the origin of 
these spirals, planets embedded in the disk may be capable of exciting the ob- 
served morphology. Assuming that this is the case, we can make predictions on 
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the locations and, possibly, the masses of the unseen planets. Such planets may 
be detected bv future multi- wavelengths observations. 

Subject headings: circumstellar matter instrumentation: high angular reso- 

lution polarization protoplanetary disks stars: individual (SAC) 206462, 
1 1 D 1 35344 B } waves 


1. Introduction 


Dynamical processes in protoplanetary disks such as turbulence or disk-planet interac- 
tion are important in understanding physical condition and evolution of disks, and planet for- 
mation processes. High resolution, direct imaging observations of cireurnstellar/protoplanetaiy 
disks can reveal non-axisymmetric structures, providing insight into such dynamical processes 
(e.g., Hashimot.o et ah 201 1). 

Recent observations have identified a class of protoplanetary disks harboring tens of 
AU-scale holes/gaps at their centers: so-called “transitional" disks. One well-studied system 
of that class is the rapidly rotating Her big F star, SAG 206462 (HD 1 35344 B. F4Ve, d = 
142 ± 27 pc, M = Muller et al 2011), The observations of CO line profiles 

(Dent et al. 2005; Pontoppidaii et al, 2008; Lyo et al. 2011) and stellar rotation (Muller 
et al, 2011) consistently indicate an almost face-on geometry ( i ™ 11°). The gap in the 
disk was predicted from I he infrared ( I R ) spectral energy distribution (SED. Brown et al. 
2007), and was subsequent ly imaged in sub- mm dust continuum at ~ Q?5 x (>725 resolution 
(Brown et al. 2009), Andrews et al. (2011) estimate the gap radius to be ~ 46 AU and the 
surface density within the gap to be 10~ fj ' 1 2 times smaller than that extrapolated from the 
outer disk. The gas in Keplerian motion surrounding the gap region is also imaged by CO 
lines {Lyo et al. 2011). The CO rovi brat ion al line observations ( Pont oppidan et al. 2008) 
and [OI] spectral line observations (van der Idas et al. 2008) indicate the presence of a gas 
disk in the vicinity (several AU-scale) of the star. SED modeling (Grady el al. 2009) and 
NER i nt e rf e i x >i net. ry ( Fede 1 e et a 1.21 )08 ) i nd i ( ■ at e 1 1 le p i ese nee of an inner dust belt, w h i cii 
is temporally variable (Silko et al. 2011] and not copianar with the outer disk (Benisty 
2011. private communication). New imaging with high spatial resolution and sensitivity is 
required to understand the inner structures of the disk. The outer portions of gaps can now 
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be resolved using 8-10 m ground-based telescopes at near infrared (NIR) wavelengths (e,g,. 
Thahnann et al 2010, for LkCa 15). 

In this Letter, we present //-band polarized intensity [PI) observations of the disk of 
SAO 206462 down to r ~ C//2 (~ 28 AU) scale at 0?06 (~ 8.4 AU) resolution. Interior to 
the sub-min resolved gap, we find spiral structures, indicative of dynamical processes. We 
use the spiral density wave theory to interpret the structure, and estimate disk's physical 
parameters. 


2, Observations & Data Reduction 
2,1. HiCIAO Observations 

SAO 206162 was observed in the H - band (1.6 pm) using the high-contrast imaging 
instrument HiCIAO (Tauiura et al. 2006; Hodapp et aL 2008; Suzuki et, al. 2010) on the 
Subaru Telescope on 2011 May 20 UT as par! of Strategic Explorations of Exoplanets and 
Disks with Subaru (SEEDS, Tamura 2009), The adaptive optics system (A0188; Hayano et 
al, 2004; Minowa et al 2010) provided a stable stellar point spread function (PSF. FWHM = 
We used a combined angular differential imaging (ADI) and polarizat ion differential 
imaging (PDI) mode with a field of view of 10* x 20* and a pixel scale of 9,5 mas pixel -1 . A 
(//3-diameter circular occulting mask was used to suppress the bright stellar halo. The half- 
wave plates were placed to four angular positions from 0°, 45°, 22.5°, and 67, 5 C in sequence 
with one 30 sec exposure per wave plate position. The total integration time of the polarized 
intensity (PI) image was 780 sec after removing low quality images with large FWHMs by 
careful inspections of the stellar PSF. 


2,2. PDI Data Reduction 

The raw images were corrected using IRAK 2 for dark current and flat-field following 
the standard reduction scheme. We applied a distortion correction using globular cluster 
M5 data taken within a few days, using IRAK packages GEOMAP and GEQTRAM* Stokes 
(QM) parameters and the PI image were obtained in t he standard approach (e.g.. Hinkley 
et al. 2009) as follows. By subtracting two images of extraordinary- and ordinary-rays at 
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each wave plate position, we obtained +Q, -Q* +! \ and — V images, from which 2 Q and 2U 
images were made by another subtraction to eliminate remaining aberration. PI was then 
given by PI = sjQ 1 + C 2 , Instrumental polarization of HiCIAO at the Nasmyth instrument 
was corrected by following .loos et ah (2008). From frame-by-frame deviations, the typical 
error of surface brightness (SB) was estimated to be ^ 5% at r ^ 0T5 when averaged over 
5 x 5 pixels (** PSF scale). Comparing different data reduc tion methods (frame selections 
and instrumental polarization estimates), we expect that the systematic uncertainty of the 
SB of PI to be - 10%, 


2.3, Contemporaneous Photometry 

Since SAC) 206462 shows variability In NIR wavelengths (Sitko et al, 2011), it is impor- 
tant to take photometry simultaneously with disk observations, //-band photometry was 
obtained just before and after the disk imaging without the coronagraphic spot with the 
adaptive optics, by sixteen 1,5 sec exposures at four spatially dithered positions. An ND10 
filter (9,8 ± 0.1% transmission) was used to avoid saturation. Using the MKO filter set. 
the H - band (A rf f 1.615 pm, FWHM = 0.29 pm, Tokunaga et al. 2002) magnitude was 
6.96 ± 0.07 mag. 

Broad-band VRIJHK photometry was obtained on 2011 May 23-26, starting within 48 
hours of the HiCIAO observation, using the Rapid Eye Mount (REM) Telescope at La Silla, 
Chile (Covino et al. 2004). The REM //- band filter has A e a — 1.65 pm. FWHM = 0.35 pm: 
broader and displaced to longer wavelengths than the MKO filter. The observed data were 
reduced differentially using SAG 206463 (AOV). The IR excess due to the inner disk (Figure 
1) was average for the range observed in 2009-201 1 (Sitko et al. 201 1 ). No significant variation 
was observed during the 2011 May observations, except for the small long-term fading trend 
(<Sin = 0.08 i 0.02 mag) over the observation period. 

Figure 1 also displays spectra obtained with the SpeX spectrograph (Rayner et al. 2009) 
on NASA's Infrared Telescope Facility {I R I F). The observations were obtained in the cross- 
dispersed (XD) echelle mode between 0,8 - 5,1 pm using a 0?8 slit (R ~ 900) and calibrated 
using HD 1 29685 ( A0\" ) with SpeXtooI (Vncca et al, 2003: Cushing et al. 2004). The absolute 
flux calibration, to correct for light loss at the spectrograph slits, was accomplished in two 
ways: using photometry and wide-slit spectroscopy (see Sitko et al, 2tJl 1), The March SpeX 
data were normalized using the REM photometry, obtained in ihe flays immediately after 
the SpeX observations, and when the star was photometrically stable. In July. SAC) 206462 
and the calibration star were observed with the SpeX prism using a 3"0 slit, which, under 
good seeing and transparency conditions, produces absolute flux to ^ 5% accuracy. The 
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REM photometry at //-band in May is ~ 0.25 mag brighter than the Subaru data. The 
Subaru photometry shows a low value even when considering the filter difference, suggesting 
that the outer disk is illuminated efficiently. 


3. Results 
3-1* Spiral Structure 

The SAO 206462 disk can be traced in PI from 072 to ^ 170 (28 — 1 10 A l ). similar to 
the range imaged with ii ST /N I CM OS (Grady et ak 2009), but with a factor of — 4 greater 
angular resolution. The total PI is 9.87 mJy ± 0.06% at 072 < r < 170. which is 0.6% of 
the stellar intensity. The total Pf at 0742 < r < 170 is 3.94 mJy ± 0.1% while the total 
intensity by HST/NICMOS F110W is 9.7 inJy (Grady et al. 2009). The average SB of PI 
at r — (7/16 is ^ 6 mJy/asec 2 , whereas the total intensity by HST/NICMOS F160W is 
30 mJy/asec 2 (Grady et al. 2009). Given the NICMOS data uncertainties, the polarization 
fraction is ^ 20 — 40% : assuming no PSF halo in the HiCIAO data and no variable self- 
shadowing/ illumination in the disk. Our measured polarization fraction is similar to that of 
HD 100546 ( 14%%g^ . Quanz et al. 2011a) and AB Aur 25% — 15%. Perrin et al. 2009). 

Figure 2 shows the PI image. The region interior to 074 is not a void and we do not see 
clear structural evidence of the cavity wall in Andrews et al. (201 1) model {R CAV = 16 AU ~ 
(y/33)* We see spiral arcs SI (east ) and S2 (south-west). The PI at the location of the spirals 
is ^ 30% larger than that extrapolated from the smoot h outer profile (bottom of Figure 2). 
The brightest port ions of the spirals roughly coincide with the bright thermal emission peaks 
at 12 /mi (Marinas et ak 2011) and lie inside the ring noted by Qoucet et ah (2006), It is 
also noted that we see a dip in PI in the north-west, probably due to the depolarization in 
the minor axis direction (see below), and that we do not see large-scale, localized shadow 
that might be cast by the inner dust belt if highly inclined relative to the outer disk. 


3,2. Azimuthal and Radial Profiles 

Here, we summarize the overall disk structure exterior to the spirals. Figure 3 (top 
panels) shows the azimuthal SB profiles at r — 075 and r — 1)77. At r ^ [7/5, SB has maxima 
around position angle (PA) ^ 50 c —60 o and 230° — 2 10°. Since the polarization is maximized 
at ^ 90° scattering (e.g„ Graham et al. 2007), it is implied that the disk major axis is at 
PA ^ 50 c -60 c . comparable to estimates by CO observations: PA - 56° ±2° by Pontoppidau 
et ak (2008) and 64° ±2° by Lyo et ak (2011), We adopt PA 55° for the major axis and 
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/ — 11 for inclination {see Section 1). Our spiral model fitting results (next Section) are 
Uti le affected even if we assume a face-on geometry. 

From CO observations, it is known that the south-west side is receding (e.g.. Lvo et al. 
201 1). Therefore, either the north-west or south-east side is the near side. We do not see an 
obvious forward scat tering excess in the NIK image. However, since the spirals are typically 
trailing, it is inferred that the south-east is t lie near side. 

Figure 3 also shows the radial PI profiles along the major axis, which is roughly consis- 
tent with in the outer part, indicating a flat (not highly flared) disk (Whitney and Hart- 
mann 1992), The radial slopes vary as PA from ^ —2 to 4.5 (fitting at, 0?6 < r < 1"G) or 

from ^ “2.5 to ~ —5 (fitting at O'.' 3 < r < 0?9) T with shallower slopes typically appearing in 
the minor axis directions. However, r“ :i is representative on average. This slope is observed 
in several other HAeBe disks {e.g., Fukagawa et ah 2010, for total intensity data), although 
HD 97048 disk exhibits a shallower slope (Quanz et al. 2011b). 


4, Spiral Structure Modeling 


Among several features in the image, the most interesting one is the aon-axjsymmetric 
spirals. In order to understand them, we propose a model based on the spiral density wave 
theory (e.g., Lin and Shu 1964: Goldreich and Tremaine 1978, 1979; Ogilvie and Lubow 
2002), assuming that NIR emission traces the disk surface density structure. With such a 
model, the spiral structures can be used to infer the disk temperature, independently of. for 
example, CO line observations. 


The shape of the spiral density wave is determined by the location of the launching 
point (corotation radius r € ) and disk's thermal and rotation profiles. When the disk rotation 
angular frequency is Q(r) oc r~ a and the sound speed is c(r) x r ' . the shape of the wave 
far from r c is given by 
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in I lie polar coordinate (r, 9\ where h v = c(r,.}/r r Q(r c ) denotes the disk aspect ratio at 
r ” /y and gives the phase. Equation (1) well approximates the shape of the density 
wave given bv the WKB theory (Rafikov 2002; Mute et al, 201 1 ), When the spiral is excited 
by a planet in a circular orbit, its location is ^ Equation (1) has five parameters. 

(r t ., tfn./v.n. J), 


9 


Two uon-axtay m metric features, SI and S2, (Figure 1) are identified ns follows. First , 
local maxima in the radial SB profiles normalized by r 2 (to take into account the dilution 
of the stellar flux) are traced at every 1° step with data at 170° < PA < 360° (SI) and at 
50° < PA < 190° (S2). The points near the minor axis (PA ^ 325°) are excluded because 
the structure there may be affected by depolarization. The points at PA > 200° (SI) and 
PA > 1 1-1° (S2) may be a part of axisyinmetrie rings since they are found to have constant 
radii. After removing these points, we have 27 (Si) and 56 (S2) points as representing 
samples of nou-axisymmetric spirals, with the opening angle of ^ ]~f for both SI and S2, 
We estimate that the uncertainty of the locution of (he maxima is given by the FWHM of 
the PSF. 

In order to fit the non-axisymnietric structures by Equation 1. we fix tv and 3 at 1.5 
(Kepler rotation) and (k b respectively, as in Lyo et al (2011), while other parameters are 
varied its (07 1 < r c < 079, 0 < 0 {i < 2tM).05 < k v < 0.25). Note that different values of 
3 yield similar results. Since it is difficult to [if SI and S2 simultaneously, they are fined 
independent ly. 

The “best-fit 1 ’ parameters are (r o 0Q,ft c ) = (0739. 204°, 0.08) for SI (reduced \ 2 — 0.52) 
and (r c ,6oJi v ) (0?9 T 353°, 0.24) for S2 (reduced \ 2 = 0.31). The spiral shapes with these 
parameters are shown in Figure 4. However, the parameter degeneracy is significant. Figure 
5 shows \ he parameter space of (r ci 0 o ) with 63.8% confidence level for h r — 0. 1 and h c — 0.2. 
Note that in Figure 5, the "best-fit"’ of (r r ,.0o) is outside the domain of confidence in some 
cases because h c is not the same as the best-fit. Despite the parameter degeneracy, the 
values of the aspect ratio which fit the shape of the spiral (h c ** 0.1) are consistent with 
those obtained from the sub-mni map of the disk (e.g.. h = 0.096(r/100AU)°‘ 15 ; Andrews et 
al. 2011). 

The spiral density wave theory predicts that the pattern speed deviates from the local 
Kepler speed; 

SW„ = 0.8 (^jj) ‘ J ' 2 (j4|:) ' [deg/yr] (2) 

is not necessarily equal to ft(r). When /, 0?5{~ 70 AL ), the spiral will move ^ 10° in a 

decade, corresponding to movement of 07 L Considering the PSF scale of our observations 
and the locations of the spirals, such deviations can be detectable over a couple of decades. 
Moreover, if the two spirals have distinct corotation radii, their relative locations change in 
time due to the pattern speed difference. Such measurements will confirm that the observed 
feature is really the density wave, providing indisputable evidence of dynamical activity. 

Note that it would be difficult to detect spirals in colder disks (smaller h e ). where spirals 
are more tightly wound, due to the blurring by the PSF. The lower detectable limit of h c 
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is typically h v ~ 0.01 - 0.03 for our set of parameters. The combination of high angular 
resolution and warm temperatures allows the spiral structure in the SAG 206462 disk to be 
resolved* Further spirals might be detectable in similarly warm disks. 


5* Summary and Discussion 

In this Letter* we present a high-resolution image of the SAC) 206462 transitional disk 
using Subaru/HiCIAO. with an inner working angle of (f!2. We discover uon-axisymmetric 
spiral features, which can be explained by the spiral density wave theory with a reasonable 
value of the disk aspect ratio (h c ™ 0,1), The model is robust in a sense that it does not 
assume the origin of such structure explicitly. 

The detection of scattered light within the sub-nun cavity itself is interesting, since 
Andrews et aL (2011) predicts tlial the sub- mm cavity is heavily depleted* Our data in 
tandem with the millimeter data may suggest that the depletion of grains at different sizes is 
not uniform. Dong et al. (2012) discuss such discrepancies between sub-mni dust continuum 
images and MR scattered light images in terms of differing spatial distributions as a function 
of grain size from a general t heoretical perspective. 

Our major assumption is that PI at FI- band, tracing the scattered light at the disk sur- 
face. actually traces the surface density variations. This assumption is valid when the disk 
is in vertical, isothermal hydrostatic equilibrium without rapid radial surface density varia- 
tions (e.g., Muto 2011). Structures near the midplane are, however, preferentially observed 
at longer wavelengths with high spatial resolution: Atacama Large Millimeter /Submillimeter 
Array (ALMA) can be an ideal instrument. 

Among several possible causes for the spiral structures (see also Hashimoto et al. 2011), 
one interesting idea is that planets excite them. In this case, the domain of possible locations 
of the planets are given by the green curves in Figure 5. If the two spirals have distinct 
corotation radii, there may be two (unseen) planets embedded in the disk. The amplitude 
of the surface density perturbation scales with the planet mass as 52/2 ^ GM p il/c z for 
tion-gap-opening low-mass planets (GM p Q/c 3 <: L e.g., Tanaka et ah 2002)* In mu data, 
the amplitude of the spiral wave is typically ^ 30% (Figure 2)* implying M p ^ 0.54/j- 

The typical error of PI in our image is ~ 5%. If this is typical of Hi Cl AO. it is capable 
of detecting the indirect signatures of planets down to M p ** 0.054/ j. ADI is promising in 
finding a point source; however, small field rotation 13°) due to the southerly declination 
(—39°) of SAG 206462 makes obtaining the total intensity difficult in our data. Future L- 
band observations may reveal thermal emission from a planet, if it exists, or its surrounding 
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(accreting) gas. 
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Fig. 1 . SED for SAO 206462 obtained by REM, adapted, in part, from Sitko et al. (2011). 
The REM observations consist of \ llIJlfK photometry. Also shown are spectra obtained 
with the SpeX spectrograph. See Section 2.3 for the data reduction techniques. 
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Fig* 2. Top: PI image of SAC) 206462 in the north- up configuration with log-stretch color 
scales* The filled orange circles at the center indicate the mask size (r = 0? 15}* The circles 
around have r — (Jf! 2. exterior to which the features are considered to he real* The right 
panel is central region's closeup. Different color scales are used to enhance the spirals labeled 
as "SD and ;i S2 r . The 'Dip' may be due to depolarization* Bottom: PI profile along the 
red line in the top right panel* The arrow indicates the location of 32. The position errors 
are uot shown for visibility. 
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Fig. 3. Top: Azimuthal SB profile at r = Of! 5 (left) and r (f!7 (right) with PA measured 
from north to east. Bottom: Radial profile along the major axis. Position errors indicate 
the FWHM of t he PSF, SB errors are estimated using f rame- \ >y-frame deviations. 
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Fig, T Left : Red and green points indicate the locations of the maxima of radial profiles 
for SJ (top) and S2 (bottom). The green ones are used for the fitting by Equation (1) 
(Section 4), Dashed lines indicate circles with r = C//39 (top) and r = Q?3 (bottom). Right: 
The spiral shape given by Equation (1) with the best-fit parameters of SI (top) and S2 
(bottom). Crosses show the locations of (r c . #q). 
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Fig. 5. — The parameter degeneracy of the fitting for SI (top) and S2 (bottom) with h c = 0.1 
(left) and A,. - 0.2 (right). The red points show a part of the data fitted. The green curves 
indicate the 63.8% confidence level for the locations of (r c . ()„). The magenta crosses are 
(r c . 0 O ) for the best- fit parameters. Note that the best -fit values of h c are different from 
those shown in the figure. 
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